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SYSTEM AND METHOD FOR SCANNING NEAR-FIELD 
OPTICAL TOMOGRAPHY 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

This invention relates to tomography and, more particularly, to near-field 
microscopy wherein an image of an object is directly reconstructed with sub-wavelength 
resolution. 

2. Description of the Background Art 

Near-field scanning optical microscopy (NSOM) is a technique to obtain 
images of surfaces with sub-wavelength resolution. The technique is particularly 
important for imaging structures where spectroscopic concerns or sample handling 
requirements dictate the use of lower frequency fields and yet high spatial resolution is 
still required. Applications range from the inspection of organic and biological samples 
to semiconductors. Various modalities are in practical use. Two prominent examples of 
such modalities are "illumination mode NSOM" and "collection mode NSOM". In 
illumination mode NSOM, a tapered fiber probe with a sub-wavelength size aperture 
serves as a source of illumination in the near-zone of the sample. The scattered field 
intensity is then measured and recorded as a function of the probe position while the 
probe is scanned over the sample. In collection mode NSOM, the fiber probe serves to 
detect the scattered field in the near-zone as the sample is illuminated by a source in the 
far zone. 
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There are certain limitations of current NSOM techniques. Despite the 
fact that the sample may present a complicated three-dimensional structure, NSOM 
produces only a two-dimensional image of the sample. Indeed, rather than being an 
imaging method, it is more accurate to say that NSOM maps the sub-wavelength 
structure of the optical near-field intensity. Under certain simplifying assumptions, such 
as homogeneity of the bulk optical properties of the sample, the maps produced in these 
experiments may be related to the sample structure. However, for the more general case 
in which the topography of the sample and the bulk optical properties both vary, the 
connection between the near-field intensity and the sample structure can be ambiguous. 

To resolve this ambiguity it is desirable to solve the inverse scattering 
problem (ISP). The ISP may be characterized as reconstructing the three-dimensional 
object structure, in this case the dielectric susceptibility of the sample, from 
measurements of the scattered field. By solving the ISP, two main issues of the prior art 
are resolved, namely, the ambiguity in connecting the sample properties and the 
measured data is resolved, and simultaneously three-dimensional, tomographic images of 
the sample are obtained. 

Historically, solving the ISP for other scattering modalities has greatly 
expanded the functionality of existing methods. For instance, pioneering analysis of X- 
ray diffraction made modern crystallography a reality. Another pioneering work brought 
medical imaging out of the era of projection radiography and into the era of computed 
tomography (CT). In any ISP the first step is to obtain a physically reasonable forward 
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model for the scattering process. For instance, in CT, a geometric model of propagation, 
neglecting any scattering, sufficiently describes the experiment. Likewise, as scattering 
becomes important, the first Born approximation is a reasonable model. By considering 
the far-zone scattered field, one may obtain a readily soluble ISP, now generally known 
as diffraction tomography (DT). The crucial step in the DT solution is to obtain a 
linearized relationship between the sample properties and the scattered field or some 
simple function of the scattered field. This may be accomplished by making use of 
models such as the first Born or first Rytov approximations. Solutions for the non- 
linear ISP may be obtained as well, but in general they suffer from mathematical 
pathologies involving convergence. 

The far-zone ISP has an inherent resolution limit imposed by the 
wavelength of the probe field. This limit may be traced to the fact that only the 
homogeneous part of the scattered field contributes to the far-zone. While in principle a 
higher resolution image may be obtained by mathematical extrapolation, this approach is 
exponentially sensitive to errors in the scattering data. However, the evanescent waves 
which contribute to the near-field carry the higher spatial frequency information about the 
scatterer. It is known that inclusion of the evanescent waves in the standard back 
propagation algorithm of two-dimensional DT enhances the resolving power of that 
method. 

The sub-wavelength resolution obtained in NSOM arises because direct 
access to the evanescent scattered waves by probing the near-field is effected. It is this 
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part of the scattered field on which the sub-wavelength structure of the scattering object 
is encoded. The NSOM ISP is thus of great interest because it obtains sub-wavelength 
resolved three-dimensional reconstructions. This technique, in accordance with the 
present invention, is referred to this as "scanning near-field optical tomography". 

Representative of the art in this technological area are the following U.S. 
Patents: (a) Fiber Optic Probe for Near Field Optical Microscopy (No. 5,485,536); (b) 
Method and Apparatus for Performing Near-Field Optical Microscopy (No. 5,479,024); 
(c) Near Field Scanning Tunneling Optical Microscope (No. 5,382,789; and (d) Near 
Field Optical Microscopic Examination of a Biological Specimen (No. 5,286,970). 

The art is devoid of any teachings or suggestions for treating the inverse 
scattering problem which is applicable to the near- field case. 

SUMMARY OF THE INVENTION 

These shortcomings, as well as other limitations and deficiencies are 
obviated, in accordance with the present invention, by devising explicit inversion 
formulas that are applicable to near- field scanning optical tomography and which provide 
a direct reconstruction of a scattering potential associated with a sample being scanned. 

In accordance with a broad method aspect of the present invention related 
to scalar waves, an image of an object is generated by: (a) probing the object with 
incident scalar waves; (b) detecting scattered waves from the object, wherein the 
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scattered waves are detected in a near-field collection mode; and (c) reconstructing the 
tomographic image by executing a prescribed mathematical algorithm with reference to 
the incident scalar waves and the scattered waves to generate the tomographic image with 
sub-wavelength spatial resolution. 

In accordance with another broad method aspect of the present invention 
related to scalar waves, an image of an object is generated by: (a) probing the object with 
incident scalar waves, wherein the incident scalar waves are generated in a near-field 
illumination mode; (b) detecting scattered waves from the object, wherein the scattered 
waves are detected in the far- field of the object; and (c) reconstructing the tomographic 
image by executing a prescribed mathematical algorithm with reference to the incident 
scalar waves and the scattered waves to generate the tomographic image with sub- 
wavelength spatial resolution. 

In accordance with yet another broad method aspect of the present 
invention related to scalar waves, an image of an object is generated by: (a) probing the 
object with incident scalar waves, wherein the incident scalar waves are generated in a 
near- field illumination mode; (b) detecting scattered waves from the object, wherein the 
scattered waves are detected in a near- field collection mode; and (c) reconstructing the 
tomographic image by executing a prescribed mathematical algorithm with reference to 
the incident scalar waves and the scattered waves to generate the image with sub- 
wavelength spatial resolution. 
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In accordance with a broad method aspect of the present invention related 
to electromagnetic waves, an image of an object is generated by: (a) probing the object 
with incident electromagnetic waves; (b) detecting scattered waves from the object, 
wherein the scattered waves are detected in a near-field collection mode; and (c) 
5 reconstructing the tomographic image by executing a prescribed mathematical algorithm 
with reference to the incident electromagnetic waves and the scattered waves to generate 
the tomographic image with sub-wavelength spatial resolution. 

In accordance with another broad method aspect of the present invention 
related to electromagnetic waves, an image of an object is generated by: (a) probing the 
object with incident electromagnetic waves, wherein the incident electromagnetic waves 
are generated in a near-field illumination mode; (b) detecting scattered waves from the 
object, wherein the scattered waves are detected in the far-field of the object; and (c) 
reconstructing the tomographic image by executing a prescribed mathematical algorithm 
with reference to the incident electromagnetic waves and the scattered waves to generate 
the tomographic image with sub-wavelength spatial resolution. 

In accordance with yet another broad method aspect of the present 
invention related to electromagnetic waves, an image of an object is generated by: (a) 
20 probing the object with incident electromagnetic waves, wherein the incident 

electromagnetic waves are generated in a near-field illumination mode; (b) detecting 
scattered waves from the object, wherein the scattered waves are detected in a near-field 
collection mode; and (c) reconstructing the tomographic image by executing a prescribed 
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mathematical algorithm with reference to the incident electromagnetic waves and the 
scattered waves to generate the image with sub-wavelength spatial resolution. 

Broad system aspects of the present invention are commensurate with 
these method aspects. 

BRIEF DESCRMTON OF THE DRAWING 

FIG. 1 is a pictorial representation for the single probe illumination mode 
showing an object under test; 

FIG. 2 is a pictorial representation for the single probe collection mode showing 
an object under test; 

FIG. 3 is a pictorial representation for the two-probe mode showing an object 
under test; 

FIG. 4 depicts a system for probing the object with waves to produce scattered 
waves which are then used to directly reconstruct the image; and 

FIG. 5 is a flow diagram for directly reconstructing the image of the object 
excited by a probing wave. 
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DETAILED DESCRIPTION 

1. OVERVIEW 

When both the phase and amplitude of the optical near- field are 
measurable, the near- field ISP may be solved, and robust inversion algorithms for both 
the scalar and vector near- field scattered waves may be obtained. In this detailed 
description, three example modalities are treated. The treatment of the vector case is 
particularly important in the near-field ISP because polarization effects are somewhat 
more complicated than for the far-field problem. The solution in all cases is derived by a 
singular value decomposition (SVD) analysis of the linearized scattering problem. The 
SVD is a generalized mode decomposition that offers considerable insight into the 
scattering problem. 

In Section 2, a discussion of near-field scattering for scalar wave-fields 
obeying the scalar reduced wave equation is presented. Specifically, expressions for the 
scattered field as measured in three distinct modalities are obtained. The first two 
modalities are variations on the illumination mode and collection mode techniques. It is 
presumed that the field may be measured, that is, the measurements are phase sensitive. 
Moreover, the third modality requires the use of two near-field probes. The three 
expressions derived for the scattered fields are shown to be special cases of a single 
expression relating the measured data to the structure of the scattering object to be 
imaged. 
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In Section 3, a unified form for all three modalities is treated and it is 
shown that, in all cases, a unique solution of the inverse scattering problem exists. 
Furthermore, the S VD for the scattering kernel is then developed and the inversion 
formula for the ISP is obtained. 

In Section 4 the vector theory of near- field scattering is presented. The 
three modalities are again discussed, as well as the unified form of the modalities. 

In Section 5, the existence and uniqueness of the solution to the inverse 
scattering problem is demonstrated and the SVD is developed. It will be seen that though 
the vector nature of light makes a significant difference in the physics of the problem, the 
underlying mathematical form is much like the scalar case and the results follow from the 
scalar case in a straight- forward manner. 

2. FUNCTION THEORETIC BASIS FOR SCALAR SCATTERING 

To understand the principles in accordance with the present invention, the 
case wherein a monochromatic scalar plane wave of frequency co incident on a medium 
characterized by a scattering potential rj(r) is first considered. The field U (r) of the 
wave satisfies the reduced wave equation 

V 2 U(r) + k 2 0 U(r) = -4^ 0 2 77(r)[/(r) (1) 
where k 2 0 = 2 n I X = co/c is the free-space wavenumber of the incident wave. The field is 
then composed of two parts, namely, 
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U(r) = t/ I (r) + ^(r), (2) 
with the incident field E/.(r) satisfying the free-space reduced wave equation, 

^U^+k^Ufi) = 0, (3) 
and the scattered field U s (r) satisfying the equation 

V 2 U s (r) + k 2 0 U 5 (r) = -4< 2 77(r)t/(r). (4) 

Following standard procedures, equation (4) may be recast as the integral equation 

C/ s (r)=^ 2 pVG(r-r')/ ? (r , )t/(r') ? (5) 

where 

itt 0 |r-r'| 

G(r-r')=f- f— (6) 
4/r j r — r f | 

is the outgoing-wave Green's function for equation (4) which obeys 

V 2 G(r) + £ 0 2 G(r) = -4*S(t). (7) 

If the scattered field is much weaker than the incident wave, the first Born approximation 
applies and the following is obtained from equation (5): 



[/,(r)=^pVG(r-r')^)t/ / (r > ). (8) 



In the discussion that follows, it will prove useful to express the Green's function G in 
the plane wave decomposition 
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G(r) = ^frf%(Qrexp[/Q-/7 + ^(Q)|z|] 5 (9) 
2;r J 

where £ Z (Q) - (k 2 0 - Q)^ 2 . The modes appearing in equation (9) labeled by Q for which 
| Q |< k 0 are the homogeneous wave modes. The modes labeled by Q for which | Q |> k Q 
are evanescent. For these modes, k z (Q) = i | k z (Q) | so that the evanescent waves decay 
5 exponentially on propagation, with a concomitant loss of high-spatial frequency 
components in the scattered field. 

2A. Illumination Mode 
Next consider the case in which the sample is illuminated by a point 
1 0 source with unit amplitude at the position r, with 

ik 0 \r-T X \ 

um= f— . do) 

|r-ii| 

This arrangement is shown pictorially in FIG. 1 wherein probe 101 is 
representative of a source which illuminates sample 105. The illumination as reflected 
1 5 from sample 105 is depicted by reflected wave 102, and the illumination as transmitted 
by sample 105 is depicted by transmitted wave 103. 

The scattered field at a point r is then given by the expression 

^(r) = ^^VG(r-r0^(r')G(r'-r 1 ). (") 
20 It is assumed that the source lies in the plane z = 0 so that the position of the source is 
r, =(p,,0) with p { being the transverse coordinate of the source. The scatterer is 
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assumed to lie only in the half-space z > 0 and to be of finite range. If the scattered field 
is measured in the z > 0 half-space, this scattered field is referred to as the "transmission 
mode". If the scattered field is measured in the z < 0 half-space, this scattered field is 
referred to as the "reflection mode". It is now supposed that the observation point r is 
5 very far from the domain of the scatterer. Then observe that for | r |»| r'| 



where k(q) = (q, sk 2 (q)) is parallel to r. Here a = 1 for the transmission mode, and 

£ = -1 for the reflection mode. It is then seen that the scattered field may be expressed in 

terms of the scattering amplitude A{p i ,q) as 



G(r-r') 



(12) 



r 



U s (r) 



(13) 



r 



where with the use of equations (9), (1 1), and (12) it is found that 



(14) 



Now define a data function ®,(q„q 2 ) by a two-dimensional Fourier 
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transform of A(p x ,q), 




(15) 



Making use of equations (9) and (14), it is readily shown that 



£,(q,) J 



20 



- 12- 



*4 
W 

0 

111 



Schotland-2 

2B, Collection Mode 
In this modality, the sample 105 is illuminated from the far zone by an 
incident plane wave and the scattered field is detected in the near-zone by means of an 
idealized point detector. This arrangement is shown pictorially in FIG. 2 wherein probe 
5 201 is representative of a detector which collects scattered waves from sample 105. 
Sample 105 is illuminated by incident wave 202 in the reflection mode, whereas sample 
105 is illuminated by incident wave 203 in the transmission mode. 

The incident wave is taken to have unit amplitude and is of the form 

If (r) = expD'q • p - isk z (q)z], (17) 
10 where, once again, e = 1 for the transmission mode and e = -l for the reflection mode. 

m 

%| The scattered field, measured in the z = 0 plane at a point with transverse coordinate p, 

? is given by the expression 

I U,(p,0;q) = k 2 0 ^Vexp[/q-pW^(q)z]G(r , -r)^(r , ) ) (18) 

which follows from equation (8). As in the illumination mode case, a data function 
15 G> 2 (q, , q 2 ) is defined by a Fourier transform of the scattered field 

<*> 2 (q„q 2 ) = ^ JdV'^U^pAq,)- (19) 

Making use of the plane wave mode representation of the Green's function of equation 
(9), it is found that 

^1^2) = 77^7 K'expO'fo, -<\2)-p + WM2)-<*Mi)]zWr)- (20) 

*,(q 2 ) 

20 It is seen that 0 2 (q,,q 2 ) = 0> 1 (-q 2 ,-q 1 ), as expected from reciprocity. 
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2C. Two Probe Mode 
Now the case wherein both the source of illumination and the detector of 
the scattered field are in the near-zone of the scatterer is considered. The source of 
illumination is taken to be a point source in the z = 0 plane with position x x = (p v 0). By 
way of notation, the case that the scattered field is measured on the z = z d plane with 
z d = 0 is referred to as reflection mode detection, and the case that the scattered field is 
measured on some plane z = z d with z d > 0 is the transmission mode case. In either 
mode, the detector position in the z = z d plane is denoted r 2 = (p 29 z d ). This arrangement 
is shown pictorially in FIG. 3 wherein probe 301 is representative of a source of 
illumination for sample 105. Probe 302 is the detector for the reflection mode, whereas 
probe 302 is the detector for the transmission mode. 

The incident field is given by equation (10). Consequently, the scattered 
field is given by 

U s (r^ 2 ) = k 2 0 fjVG(r 2 -rMr)G(r- ri ). (21) 
Taking the Fourier transform over both spatial coordinates 

°3(qi>q2) = 7T^ U 2 p x d 2 p 2 exp[K<li- A -q 2 • PiWfa&Pi^d)* ( 22 ) 
(2/r) J 

and using the plane wave decomposition of equation (9), the following obtains: 

<p 3 (qi>q 2 )=- ,^trr x K'«pfr"(qi -^yp^^O-^Mi^M^ (23) 
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3. NEAR-FIELD INVERSE SCATTERING PROBLEM: SCALAR CASE 

In the near-field inverse scattering problem (ISP), the dielectric 
susceptibility or scattering function rj(r) is reconstructed from measurements of the 
scattered field. It may be observed that equations (16), (20) and (23) are all of similar 
form. In order to describe the inverse scattering problem for all three modalities to be 
discussed, the following data function is introduced: 

fcfai . q 2 ) = w(q p q 2 ) jd \ exp {i(q, - q 2 ) p + i[k z (q, ) - sk 2 (q 2 )] >J?(r), (24) 

where 

^(q^cb) = ^1(^1^2) for the illumination mode, (25a) 

^1^2) = *2(~q 2 rqi) fo r the collection mode, (25b) 

^(q* » q 2 ) = ^3 (qp q 2 ) for the two P robe mode * ( 25c ) 

Through appropriate choice of the weighting function w(q p q 2 ) all three data functions, 
Oj , 0 2 , and 0 3 may be simply related to O , that is, 

ik 2 

w (qi > q 2 ) = — ~ for the illumination mode, (26a) 

*,(qi) 

ik 2 

w (qi » q 2 ) = — — for the c °H ect i° n mode, and (26b) 

*,(q 2 ) 

wfai » q 2 ) = — 9 f° r the P r °b e mode. (26c) 

**(qi)*z(q 2 ) 



In all cases, it is assumed that the scattered field can be measured with 
point-like detection and that the measurement yields both amplitude and phase. 
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The immediate focus is on the derivation of the singular value 
decomposition (SVD) of the operator relating O to 77 as defined by integral equation 
(24). The SVD provides considerable insight into the structure of the forward scattering 
problem. It also provides a direct technique to construct the pseudo-inverse of the 
scattering operator. The pseudo-inverse may be used to construct the unique solution to 
the inverse problem with minimum L 2 norm. In addition to the construction of the 
minimum norm solution, the SVD presents an obvious procedure by which to regularize 
the pseudo-inverse through filtering small singular values. This is a significant point for 
the near-field ISP because of the presence of evanescent waves which exhibit exponential 
decay. This exponential behavior is a reflection of the ill-posed nature of the problem. 

To commence, write equation (24) as 

<*>(q l5 q 2 )= j^(q 1? q 2 ;r)77(r) 5 (27) 
where the forward scattering operator £(q 1? q 2 ;r) is given by 

£(qi,q 2 ; r ) = ^(qpq2) ex P^ < 28 > 

It is desired to find the SVD of AT which is of the form 

K(qpq 2 ;r) = f^<(r)G CT (q p q 2 ) (29) 
where a is the singular value associated with the singular functions F a and G a . The 
singular values and singular functions satisfy the relations 

JdVA;(q p q 2 ;r)F <r (r) = oG ff (q p q 2 ), (30) 

and 
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\d\ \d 2 q 2 K\c [{A2 -x)GMv^) = ^(v). (31) 

Equations (30) and (31) are sufficient to define the singular functions and singular values. 

It is also possible to derive from these equations the identities 

K*KF 0 = <r 2 F a , (32) 
KK'G a =a 2 G c . (33) 

Here,/: is the integral operator whose kernel is A:(q„q 2 ;r). The F a and G a are thus the 

eigenfunctions with eigenvalue cr 2 of K*K and KK\ respectively. 

To obtain the singular values and singular functions, it is worthwhile to 
first consider the restricted problem wherein q 2 is fixed and obtain the SVD of 
K(q l , q 2 ; r) which for fixed q 2 is denoted K qi (q, , r). Then the following obtains: 

K q2 ( qi ,r)= \d 2 Q*¥ff\r)gmq 2 ), (34) 
Vq'^Q' (35) 

Oq =<t q^ 2 ' (36) 

where the / Q Q2 and g^ 2 are the singular functions for fixed q 2 with singular value 
a q Q 2 , labeled by the two dimensional vector Q. The eigenvalue equations then become 
<K q /Q 2= K) 2 / Q q2 , (37) 

By inspection, it is seen that 

4 2 (q 2 ) = ^(Q + q 2 -q 1 ), (39) 
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f* 2 (r) = (of )-' w (q 2 + Q, q 2 ) exp HQ • p - ift (q 2 + Q) - A] (q 2 )M . 



(40) 



(Note that an explicit expression for Oq 2 may be obtained but is not needed.) The 
/ Q q2 are orthogonal for fixed q 2 , 

(/ Q q2 |/ Q q . 2 H(Q-Q')- ( 41 ) 

The /q 2 are also orthogonal for different q 2 , that is, 

(/ Q q2 \f^) = H^ r'A/fo^jQWQ-Q'), (42) 

which defines M(q 2 ,q' 2 ; Q). 

Returning now to the main problem of finding the SVD of AT(q 1 ,q 2 ;r) 3 

make the ansatz that the G's are constructed from the g£ : 

G QQ .(q 1 ,q 2 ) = C Q( (q 2 ;Q)^ 2 (q 1 X (43) 
where the singular values and the singular functions are labeled with two-dimensional 
vectors, Q and Q'. Making use of the relation 

KK*G QQ , = Jj 2 ^ q <G QQ . (44) 
it is found that equations (33) and (38) imply that 

p^'M(q,q I ;Q)C Q .(q , ;Q) = a 2 QQ C Q ,(q;Q). (45) 

That is, C Q .(q;Q) is an eigenfunction of M(Q) labeled by Q with eigenvalue o 2 qq-. 
Since M(Q) is Hermitian, the C Q ,(q;Q) may be taken to be orthonormal. It then follows 
from equation (31) that 

F QQ .(r) = ld 2 qC Q ,(q;Q)alf*(v). (46) 
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Thus the SVD of K(q l9 q 2 ;r) is given by 

Kiq^'j) = {t/ 2 ^ 2 e^ QQ ^Q Q .(r)G QQ ,(q 1? q 2 ). (47) 

An explicit evaluation of M(Q) shows that M(q,q';Q) is divergent when 
|q|, | q'| ,|Q+q|, and | Q + q'| are all less than k 0 and also in the case of transmission mode 
when k z (Q + q) + k 2 (Q + q' ) > k z (q) + k z (q f ). Of course, this divergence is unphysical. 
Any scatterer must be of finite thickness and so, to avoid divergencies, the integration 
over r is limited to the region 0 < z < z^. Then M(Q) takes the explicit form 

M(q, q';Q) = /4/r 2 w* (q, Q + q)w(q f ,Q 4- q') 

x 1 - exp{/(<(Q + q) + k m s (q) - sk z (Q + q') - k 2 ^))z^} 
<(Q + q) + K(q) - ek z (Q + q') - * f (q') 

In practice, specification of z^ acts as a means for including some prior knowledge 
about the sample 105 in the model. 

3A. The Pseudo-inverse 
To reconstruct rj (r) from the data function $(q p q 2 ) it is necessary to first 
solve the integral equation Krj = O. Since 77 and O belong to different Hilbert spaces, 
the inverse of K as a linear operator is not defined. Instead, a substitute for the solution 
is sought which is defined to be a minimizer of || Kt] - O || . Among all such solutions, it 
is conventional to choose the one with minimum L 2 norm. This so-called pseudo inverse 
solution /7 + (r) is unique and is given by 

- 19- 
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where A: + (r;q,,q 2 ) is the pseudo-inverse of /T(q,,q 2 ;r). The pseudo-inverse may be 
obtained from the SVD of K and is given by the following 

^ + (r;q„q 2 ) = [dQd 2 Q— F w (r)G^.(q I>q2 ). (50) 

Using the explicit expressions obtained for the singular functions, K + is found to be given 
5 by 

/r(r;q p q 2 ) = \d 2 Q\d 2 Q J^V(Q + q 2 -q,) 

x < Q .C Q .(q; Q)C;(q 2 ;Q)w*(q + Q,q) 

M x exp {-/Q • p - i(k\ (q + Q) - ek] (q))z} . (51) 

Q 

IS It may be noted that 

K 

J 10 p 2 g-^C Q .(q;Q)C;.(q 2 ;Q) = M l (q,q';Q), (52) 

^QQ' 

Sip: 

p where M _1 (q,q';Q) is the qq' matrix element of M (Q). Thus 

Q 

I // + (r)=fA4pW(fi+?2-?i)^ , (q ) ^^ 

S x exp HQ • p - i(kl(q + Q) - <(q)>}d>(q 1? q 2 ), (53) 

which is the inversion formula for the near-field ISP. 

15 



3B. Regularization 
It is known in the art that the singular values for K may be very small and 
20 the ISP is ill-posed. Accordingly, the kernel K + is highly singular and must be 
regularized to obtain a stable image reconstruction algorithm. Any of the standard 
approaches, such as Wiener filtering or the Tikhonov techniques, may be used. For the 
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SVD-derived pseudo-inverse, these methods result in modification of the singular values 
via the introduction of a suitable regularizer R(cr). More explicitly, equation (50) is can 
be written as 



^ + (r;q P q 2 ) = U 2 Qd 2 Q '— i?(<r QQ .)F QQ ,(r)G^( qi ,q 2 ). (54) 



The simplest choice of R(a) consists of simple truncation, which leads to 

for some a mn . As a variant of this approach, the regularizer may be 

#Oq Q 0 = ^qq-0(c-qq' " P max Q (cr QQ ,)), (56) 
where f3 is the regularization parameter and max e {cr QQ ,) is the largest singular value for 
a fixed Q. This has the effect of conditioning each of the matrices M(Q) individually. If 
Tikhonov regularization is used, 



/L + cTqq, 



where X is a regularization parameter. This choice leads to smoothing of 7] + (r) by 
penalizing functions with large L 2 norm. 
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4. VECTOR SCATTERING THEORY 

Now the attention is focused on the vector theory of light. The 
mathematical treatment of the vector near-field ISP is quite similar to the scalar case; 
reference is made to the scalar case as needed. Nevertheless, the full vector theory is 
necessary to consider in depth because the scalar approximation of the scattering of 
electromagnetic waves breaks down for the case of the dielectric susceptibility 7j(r) 
varying on subwavelength scales. In addition, the vector approach allows for the 
incorporation of effects due to polarization. 

The attention is on non-magnetic materials, so the focus is only on the 
electric field E. Consider a monochromatic field of frequency &. The electric field then 
satisfies the reduced wave equation 

V x V x E(r) - £ 2 E(r) = 4^: 0 2 ^(r)E(r), (5 8) 
where ^(r) is the dielectric susceptibility. The field is taken to be composed of two parts, 

E(r) = E'(r) + E*(r), (59) 
where E ! (r) is the incident field satisfying the homogeneous equation 

VxVxE^-^E^O, (60) 
and E 5 (r) is the scattered field satisfying the equation 

V x V x E s (r) - £ 2 E s (r) = 4^ 0 2 ^(r)E(r). (61) 
Equation (61) may be recast in integral form 

E«(r)-^ 2 \d>r<G4(T-t)E p VM?) 9 (62) 
where G a/3 is the electromagnetic free-space Green's tensor satisfying the equation 
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V x V x G afi (r) - k 2 0 G a , (r) = 4*J(r)<^ . (63) 
The following notation is used throughout, namely, Greek subscripts indicate vector 
components and the summation convention applies to repeated indices. 

5 Within the accuracy of the first Born approximation, the scattered field is 

given by 

Ei(r) = k] JrfVG^(r-r')4(0^). (64) 
M, The Green's tensor may be obtained from the scalar Green's function by observing that 

1 G a/} (r) = (S a/}+ k-%d,)G(r). (65) 

in 

f II 1 0 Thus it becomes possible to obtain the plane wave decomposition of G ap directly by the 

M 

expression 

G a Ar) = \d 2 QK(Q)~%^) «PflQ • P + *,(Q) \ z \)> ( 66 > 
where 7^(Q) accounts for the polarization of the modes and is given by the expression 

T a ,(Q) = S a ,-K%(Q)k,(Ql (67) 

15 and *(Q) = (Q, £*,(<»)■ 

4A. Illumination Mode 
Now turn to the modality considered in Section 2.A in which sample 105 
is illuminated by a point source with position r x = (p 1? 0) which is scanned in the plane 
20 z = 0. The incident filed is given by 



W 
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where p is the dipole moment of the source of the incident field. Using equation (64), the 
scattered field is seen to be of the form 

E; (r) = K J<* V G ap (r - r' ) V (f)G^ r (r'-r { )p r . (69) 
In the far zone of the scatterer, for | r |»| r'|, the asymptotic form of the Green's function 
is given by 

ikr 

G aP (r-n~T a ,(q)^-e-^\ (70) 

where equation (12) has been deployed and k(q) = (q, sk 2 (q)) is in the direction of r. It is 
then seen that the scattered field may be expressed in terms of the scattering amplitude 

e ikr 

E:(r)~klA ap {p lA )p p --. (71) 

Making use of equations (66), (69), and (70), it is found that 

4*(A.q) = \d^G r ^-x)tx V [-iq ■ p'-i£k z (q)z'Mn (72) 

As before, a data function is defined through a Fourier transform of the angular part of 
the far-zone scattered field, 

^M^) = ^\d 2 p^A^{ M2 )p p . (73) 
Making use of the plane-wave decomposition equation (66), it is found that 

<^Wq2)=7^Mq2)^q,)/v^ 

(74) 

It may be noted that ®^ is identical to O, except for the polarization factors, 
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«£ ) (qi.q a ) = r ((r (q a )r rf (q I )/i / ,* I (q I ,q a ). (75) 

4B. Collection Mode 
Now consider the modality discussed in Section 2B in which an incident 
5 plan wave illuminates sample 105 and the scattered field is measured in the z = 0 plane. 
The incident field may be expressed as 

K (r) = e a exp[/q • p - isk z (q)z], (76) 

|* where e a is the polarization incident field. The scattered field, measured in the z = 0 

Q 

*g plane at the position r = (p d ,0) , is given by the expression 

iff, 

| 10 ^(pAq) = k\s p J ^ /VG a/? (r-r•)exp[ I ■q•pW^ ^ (q)z•]^(r , ). (77) 

!^ Now define the data function as the Fourier transform of the scattered field 

<&2(q,.q 2 ) =^ J^V iq2 '^(PAq,), (78) 
and use the plane wave decomposition equation (66) to obtain the explicit form 



s 

p. 



J. 

-or^'-'" 27V 



m 
0 



&?{^) = ^T ap ^ 2 )e p \d'rsM^^ (79) 

1 5 Again it may be seen that cD^ 2) and 0 2 are simply related by the expression 

^ 2) (qpq2) = ^(q 2 )^0 2 (q 1 ,q 2 ). (80) 

4C. Two Probe Mode 
Now the case wherein both the source of illumination and the detector of 
20 the scattered field are in the near-zone of the scatterer is considered. The source of 
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illumination is taken to be a point source in the z = 0 plane with position r, = (p x ,0). By 
way of notation, the case that the scattered field measured on the z = 0 plane is referred 
to as the reflection mode, and the case that the scattered field measured on the 
plane z = z d with z d > 0 is referred to as the transmission mode case. In either mode, the 
detector position in the z = z d plane is denoted r 2 = (p 2 ,z d ). The incident field is given 
by equation (68). Consequently, the scattered field is given by 

K(Pv®>Pi> z d) = K ftr'G^ -r')G^(r'-r 1 )^(r*)- (81) 
Taking the Fourier transform over both spatial coordinates defines a data function 

®a\w 2 )=7j-J f</W 2 exp[/(q, - a -q 2 ■ p 2 )}E s a {p x ,0; p 2 ,z d ). (82) 

(27C) 

It is then found that the following obtains: 

i 4 i£k z {q 2 )zd 

x f</Vexp{*(q, -q 2 )- p + i[kMi)-^A<l2)] z Mr)- (83) 

The data functions <D^ 3> and 0> 3 are simply related by 

^l 3) (qpq2) = 7; /? (q 2 )^(q 1 )^03(q 1 ,q 2 ). (84) 

5. NEAR-FIELD INVERSE SCATTERING PROBLEM: VECTOR CASE 

As in the scalar case, the inverse problem is treated for all three modalities 
by considering a data function of the form 

O a (q„q 2 ) = w a (q„q 2 ) J</V exp^q, -q 2 ) • p + /[^(q,) -^ z (q 2 )]z}/7(r), (85) 

where 



-26- 



Schotland-2 

< £ a (<lp<l2) = *a ) ( ( li> £ h)> for illumination mode (86a) 
^3(^1^2) = *a > (- < l2'- c li)' for collection mode (86b) 
® a (. ( li,<l2) = $a ) ( c lv < ±2)> for two-probe mode. (86c) 

Through appropriate choice of the weighting function w a (q p q 2 ) all three data functions 

may be simply related to O a , that is, 

wA^Q 2 ) = --^-JaM2) T A^Pp > illumination mode (87a) 
pfc) = -T^T^iq^ , collection mode (87b) 

w AVi><l2) = 7T, T ap {q 2 )T pr {^)p r , two-probe mode. (87c) 

The vector integral equation (85) differs from its scalar counterpart of equation (24) only 
by a factor associated with polarization. It is clear then that by measuring a fixed 
component of the scattered field that the scalar inversion formula of equation (49) may be 
used to reconstruct ^(r). It is advantageous, however, to utilize measurements of all 
three components of the scattered field. In this case, the ISP is over-determined and the 
SVDofthe operator relating O a and 77 is needed. 

To obtain the required SVD, begin by writing equation (85) as 

0 o (q 1 ,q 2 )=|^ a (q 1 ,q 2 ;rMr), (88) 

where the forward scattering operator K a (q„q 2 ;r) for the vector wave case is given by 
K a (q, , q 2 ; r) = w a (q, , q 2 ) exp {i(q, - q 2 ) • p + i[k z (q, ) - sk z (q 2 )]z} . (89) 
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Consider first the restricted problem where q 2 and a are fixed, and obtain 
the SVD of the vector operator with components ^ a (q„q 2 ;r) which, for fixed q 2 and a , 
is denoted K q2 (q„r). The SVD is then obtained as follows: 

* q2 *(q.,r)= iSQo&fgtog^), (90) 
where the f$ 2 are the singular functions for fixed q 2 and a with singular values cr% a . It 
is readily seen that 

/q« =(<)"'w (q 2 +Q,q 2 )exp{-/Q-p-p;(q 2 +Q)-^(q 2 )M- (91) 
As before, an explicit expression for <r% a is not required. The f^ a obey the 
orthogonality conditions 

(/ q *|/q4W(Q-Q')' (92) 

and 

{fo 2 a | f$i ) = (« r^aM^ ; Q)*Q - Q')' (93) 

which defines M ap (q 2 , q' 2 ; Q). 

To find the SVD of the vector operator £(q p q 2 ;r), it is necessary to 
determine the eigenfunctions of KK* which satisfy 

X^^Gj 3 .(q 1 ,q 2 ) = (Tj Q .G e a e .(q I ,q 2 ), (94) 

p 

where the eigenfunctions G QQ ,(q„q 2 ) and the singular values cr QQ , have been labeled 
with two-dimensional vectors Q and Q' . Make the ansatz that the G 's are constructed 
from the : 
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G Q a Q .(qpq 2 )=Q(q2;Q)4 2 (q 1 ). ( 95 > 

Making use of equation (90), it is found that the Cg.(q 2 ;Q) are eigenfunctions of 
A/^(q,q';Q), 

jrf 2 9 'M a/? (q,q 1 ;Q)Cj(q I ;Q) = ^ Q .q.(q;Q), (96) 

5 with eigenvalues a^,. Next, note that it follows from equation (93) that 

M a/3 (q^;Q) = w a (q + Q,q)w^+Q,q')A(q^;Q), (97) 

where 

2 l-exp{/(fifc;(Q + q) + ^(q)-^(Q + q')-^(q'))^J (98) 
i4(q ' q;Q) " ^(Q + q) + *;(q)-^(Q + q')-W) 

Since the prefactor w a (q + Q,q)w /; (q , +Q,q') of ^(q,q*;Q) has the form of a projection 
1 0 operator, it is found that C£ (q; Q) may be expressed as 

w g (q+Q>q) (99) 



1 C?(«Q)-«q.Q) fl(q>Q) 



where 



. V' 2 



Q(q,Q) = 



I|w a (q + Q,q)| 2 • (100) 

V a ) 



The 0 Q .(q;Q) are easily seen to be the eigenfunctions of 
1 5 A(q, q' ; Q) = Q(q* , QM(q, q* ; Q)Q(q, Q) with eigenvalues cr 2 e . , that is, 

}j 2 ? M(q,q , ;Q)^(q';Q) = ^ Q ,(q;Q). (101) 
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A second class of eigenfunctions of M^(q,q';Q) may also be seen to exist. However, 
these eigenfunctions have vanishing singular values and do not contribute to the SVD of 
K. 

To complete the development of the SVD, the singular values of F QQ ,(r) are 
required, which as follows from equation (31), are given by 

F QQ .(r) = < 7-.2<G« Q . (102) 

a 

=<.^ 2 ?Sq,(q;QX/4(r). 

a 

Using this result, it is seen that the SVD of -£ a (q„q 2 ;r) is given by 

KMr>W)= p 2 ^ 2 eV QQ .F^.(r)G^( q] ,q 2 ). (103) 

The pseudo-inverse of K a (q, , q 2 ; r) may be obtained from the SVD of K a 
and is given be the expression 

C(r;q„q 2 )= \d 2 Qd 2 Q— F QQ .(r)G«;,(q 1 ,q 2 ). (104) 



Using the expressions previously obtained for the singular functions, K + a is found to be 
of the form 

K(x- Av q 2 )=\d 2 Q\d 2 Q JrfV(Q+q 2 -q.) 

x«(q 2 ;Q)^.(q;Q)w;(q + Q,q) 

x exp HQ-p - i(k]{q + Q) - ^(q))z} . (105) 

Using the fact that 
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w« (q + Q> (q'+Q> q') 
J^g«(q;Q)^.(q';Q)^-'(q.q';Q) ° W fl( ^ q . Q) ■ (106 ) 



and the expression for the pseudo-inverse solution 

// + (r) = Jrf 2 ?^ 2 ? 2 C(r;qpq 2 )^(qpq 2 X a 07 ) 

the following obtains: 

77 + (r)= frfV 2 ft K^ 2 ^7?r^^(Q + q2 -q 1 M" 1 (q 25 q , ;QK(q 2 + Q,<7 2 ) 

x expHQ ■ /> - i(*;(q + Q) - ^faM^i^X (108) 
which is the inversion formula for the near-field ISP. 

To this point it has been assumed that the scattered field may be measured 
in its entirety. It may the case that only certain components of the scattered field are 
measured. Let | be a unit vector in the direction in which the scattered field is 
measured. The measured quantity will then not be E s a but | • E\ The appropriate data 

function to consider for all three modalities in then 0(q l9 q 2 ) = £® a (q v q 2 )' This 
modifies the weight function without otherwise changing the analysis. 

6A. System 

As depicted in high-level block diagram form in FIG. 4, system 400 is a 
tomography system for generating an image of an object using measurements of scattered 
waves emanating from an object in response to waves illuminating the object. In 
particular, object 105 is shown as being under investigation. System 400 is composed of: 
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wave source 420 for probing the object 105; data acquisition detector 430 for detecting 
the scattering data corresponding to the scattered waves from object 105 at one or more 
locations proximate to object 105; position controller 440 for controlling the locations of 
detectors 430 and sources 420; and computer processor 450, having associated input 
5 device 460 (e.g., a keyboard) and output device 470 (e.g., a graphical display terminal). 
Computer processor 450 has as its inputs positional information from controller 440 and 
the measured scattering data from detector 430. 

Computer 450 stores a computer program which implements the direct 
reconstruction algorithm; in particular, the stored program processes the measured 
scattering data to produce the image of the object or sample under study using a 
prescribed mathematical algorithm. The algorithm is, generally, determined with 
reference to an integral operator relating the scattering data to the forward scattering 
operator (e.g., integral equation (27) wherein the forward scattering operator is 
AT(q„q 2 ;r) or by integral equation (88) wherein the forward scattering operator is 

K a (q„q 2 ;r)). 

6B. FLOW DIAGRAM 

The methodology carried out by the present invention is set forth in high- 
20 level flow diagram 500 of FIG. 5 in terms of the illustrative system embodiment shown in 
FIG. 4. With reference to FIG. 5, the processing effected by block 510 enables source 
420 and data acquisition detector 430 so as to measure the scattering data emanating from 
object 105 due to illuminating waves from source 420. These measurements are passed 
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to computer processor 450 from data acquisition detector 430 via bus 431. Next, 
processing block 520 is invoked to retrieve the pre-computed and stored forward 
scattering operator. In turn, processing block 530 is operated to execute the 
reconstruction algorithm set forth with respect to the foregoing equations, thereby 
determining the scattering potential rj{r) for the scalar case or dielectric susceptibility 
r/(r) for the electromagnetic case. Finally, as depicted by processing block 540, the 
reconstructed tomographic image corresponding to rj{r) is provided to output device 470 
in a form determined by the user; device 470 may be, for example, a display monitor or a 
more sophisticated three-dimensional display device. 

Although the present invention have been shown and described in detail 
herein, those skilled in the art can readily devise many other varied embodiments that still 
incorporate these teachings. Thus, the previous description merely illustrates the 
principles of the invention. It will thus be appreciated that those with ordinary skill in the 
art will be able to devise various arrangements which, although not explicitly described 
or shown herein, embody principles of the invention and are included within its spirit and 
scope. Furthermore, all examples and conditional language recited herein are principally 
intended expressly to be only for pedagogical purposes to aid the reader in understanding 
the principles of the invention and the concepts contributed by the inventor to furthering 
the art, and are to be construed as being without limitation to such specifically recited 
examples and conditions. Moreover, all statements herein reciting principles, aspects, 
and embodiments of the invention, as well as specific examples thereof, are intended to 
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encompass both structural and functional equivalents thereof. Additionally, it is intended 
that such equivalents include both currently know equivalents as well as equivalents 
developed in the future, that is, any elements developed that perform the function, 
regardless of structure. 

In addition, it will be appreciated by those with ordinary skill in the art that 
the block diagrams herein represent conceptual views of illustrative circuitry embodying 
the principles of the invention. 
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